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ABSTRACT 

Context. The observations of the proplyds in the Orion Nebula Cluster, showing clear evidence of ongoing photoevaporation, have provided 
a clear proof about the role of the externally induced photoevaporation in the evolution of circumstellar disks. NGC 6611 is an open cluster 
suitable to study disk photoevaporation, thanks to its large population of massive members and of stars with disk. In a previous work, we 
obtained evidence of the influence of the strong UV field generated by the massive cluster members on the evolution of disks around low-mass 
Pre-Main Sequence members. That work was based on a multi-band BVIJHK and X-ray catalog purposely compiled to select the cluster 
members with and without disk. 

Aims. In this paper we complete the list of candidate cluster members, using data at longer wavelengths obtained with Spitzer/IRAC, and we 
revisit the issue of the effects of UV radiation on the evolution of disks in NGC 661 1. 

Methods. We select the candidate members with disks of NGC 6611, in a field of view of 33' x 34' centered on the cluster, using IRAC 
color-color diagrams and suitable reddening-free color indices. Besides, using the X-ray data to select Class III cluster members, we estimate 
the disks frequency vs. the intensity of the incident radiation emitted by massive members. 

Results. We identify 458 candidate members with circumstellar disks, among which 146 had not been revealed in our previous work. 
Comparing of the various color indices we used to select the cluster members with disk, we claim that they detect the excesses due to the 
emission of the same physical region of the disk: the inner rim at the dust sublimation radius. Our new results confirm that UV radiation from 
massive stars affects the evolution of nearby circumstellar disks. 



1. Introduction 

Circumstellar disks around young Pre-Main Sequence (PMS) 
stars have been deeply studied since the d iscovery of their in- 
frared emission (Mendoza Vl ll966l. l 1968). Until now, several 
disk models have been developed and a large number of star- 
forming regions have been observed, in order to understand 
the physical propert ies and the evolutio n of circumstellar disks. 
Several authors (e.g. lHaisch et al.l2 001) claim that the emission 
from the inner region of the disks, responsible for NIR excess 
in T-Tauri stars, declines on a timescale between ~1 and ~10 
Myr. Important clues about the typical timescale of disk evolu- 
tion come from the studies of PMS stars that have cleared the 
inner region of their circumstellar disks, thus showing excesses 
in mid- and far-infrared but not in near-infrared bands. These 
stars with disk are usually considered to be in a transitional 
phase between Class II and disk-less Class III PMS stars. 



The evolution of circumstellar disks can be influenced by 
nearby massive stars, as it was shown by the observations of 
the evaporating protoplanetary disks in the Trapezium, near the 
massi ve star ©' Ori (see, for example. IStorzer & Hollenbachl 
19991 and references therein). In this case, the photoevaporation 
process, induced by UV radiation emitted by ©' Ori, dissipates 
the nearby disks on timescales shorter tha n 10 6 years. On the 
other hand, several authors (for example Eisner & Carpenter , 
2006) claim that the externally induced photoevaporation is not 
an efficient mechanism for the truncation of circumstellar disks 
around low mass PMS stars. The debate about this topic is still 
open. 

Comparing the spatial variations of disk frequencies respect 
to the positions of massive stars in young open clusters it is 
possible to study how the radiation from massive stars affects 
the evolution of circumstellar disks and the formation of new 



genera tions of stars. Bal og et al.l (2007) and iGuarcello et al 
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(2007), hereafter GPM07, applied this approach to two young 
clusters: NGC 2244 and NGC 6611 respectively. In their 
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Spitzer/IRAC-MIPS study of NGC 2244 (distant 1.4-1.7 Kpc 



from t he Sun and with an age of about 2-3 Myr), Bal og et al 



d2007h found that disk frequency drops quickly at a distance of 
0.5 parsec from the massive cluster members. For larger dis- 
tances, however, stars with disk do not experience effects in- 
duced by massive stars, since the disk frequency is not spatially 
correlated with their positions. 

Also the young open cluster NGC 6611, in the Eagle 
Nebula, is an ideal target for this kind of study, thanks to its rich 
PMS population and its large number of massive members (56 



with spectral classes earlier than B5, Hillenbra nd et aU 11993) 



that are distributed irregularly in the central region of the clus- 
ter. In order to obtain a reliable list of cluster members GPM07 
used a membership criterion based on infrared excesses (to se- 
lect members with disk) and X-ray emission (to select mem- 
bers without disk). In fact, thanks to the high X-ray luminos- 
ity of PMs stars, this criterion allows to obtain an unbiased 
sample of PMS cluster members. With this method, a total of 
1122 candidate PMS stars associated to NGC 6611 or to the 
outer part of the Eagle Nebula were identified. In NGC 6611, 
GPM07 found evidences that the UV radiation from massive 
stars have dissipated the disks around nearby PMS members on 
short timescales, since the disks frequency declines for small 
distances from massive members. To obtain this result, GPM07 
calculated the flux emitted by all the massive stars incident on 
each cluster members and then obtained the disk frequency for 
various bins of incident flux. To this calculation, the projected 
distances were used, instead the real ones. However, GPM07 
showed that this approximation have not affected their result. 

The effects of the energetic radiation from massive stars 
on the evolution of nearby circumstellar disks are also studied 
with simulations of the evolution of open clusters. For example 
Adam s et al d2006h claim that externally induced photoevap- 
oration does no affect the evolution of disks in clusters with 
less than 1000 members. Ins t ead, m ore populated clusters, as 
shown by iFatuzzo & Adamsl (12008b . can be environments in 
which disks are efficiently evaporated in small timescales, since 
they have intense UV fields practically independent from the 
number of stellar members. 

This paper is the follow up of the GPM07 study, enrich- 
ing the list of cluster members by using Spitzer/IRAC data and 
confirming their principal results; both works, therefore, will 
be the base of a subsequent paper, focused on the analysis of 
the Spectral Energy Distributions of selected cluster members. 

The Spitzer/IRAC observations of NGC 6611 used in 
this work are obtained from the Galactic Legacy Infrared 



Mid-P lane Survey Extraordinaire (GLIMPSE; iBeniamin et al 



2003). The Infrared Ar ray Camera (IRA Q, mounted on 



Spitzer Space Telescope (IFazio et all 120041) . allows studying 



stars with circumstellar disk in four infrared bands, centered 
on 3.6 ym, 4.5 ym, 5.8 ym and 8.0 ym. In this region of the 
spectrum, the contribution from the stellar photosphere is usu- 
ally small compared to the emission from disk and envelope, 
and the effects of extinction by interstellar medium is smaller 
than at shorter wavelengths. These facts permit a more reliably 
selection of stars with circumstellar disks than using 2MASS 
data alone. 



Our paper is organized as follow: in Sect. [2] we show the 
results of the cross identification between GLIMPSE catalog 
and that compiled in GPM07; in Sect.[3]and|4]we will describe 
the identification of stars with disk respectively using IRAC 
color-color diagram and suitable reddening-free color indices 
that combine optical and NIR colors; in Sect. [5] we will com- 
pare the two diagnostics used in this paper; in Sect. [6] we de- 
scribe the catalog of the cluster members; in Sect. [7] we will 
re-examine the results obtained in GPM07 using the updated 
member list. 

2. Multi-wavelength catalog 

2.1. Description of GPM07 multiband catalog and 
NGC 6611 parameters 

GPM07 compiled a multi band catalog of NGC 66 1 1 , using op- 
tical observations in B VI bands in a 33' x 34' field of view cen- 
tered on the cluster, obtained with ESO 2.2m/WFI camera (PI: 
Mundt), 2MASS public data of the same sky region and an X- 
ray observation per formed with Chand ra/ ACIS -I in a 17' x 17' 
region (obs.ID 978: lLinskv e"taill2007h . Hereafter, we will call 
WFI FOV and ACIS FOV the fields of view of the respective 
instruments; the latter FOV is contained in the former and it is 
approximately at its center. The catalog compiled in GPM07 
consists of 38995 sources falling in the WFI FOV. 

Using this catalog, GPM07 estimated the age of the PMS 
members (that are mostly younger than 3 Myr), the distance of 
the cluster (~ 1750 parsec), the anomalous reddening law (with 
Rv ~ 3.3), the average extinction for cluster members (Ay = 
2.6"') and the relaxation time for the core (4.2 Myr, longer than 
the age of PMS members). The age of the cluster, however, is 
not well constrained as other cluster parameters. For example, 
several authors claimed th at star formation in this nebula st arted 
about 6 Myears ago (see Gvaramadze & Bomansl 2008 for a 
detailed discussion about the age spread in this cluster). 

In the ACIS FOV we found 997 X-ray sources us- 
ing PWDetect, a wav elet-based source detection algorithm 



(Da miani et al.L 11997b . with 10 expected spurious sources. 



Among the X-ray sources with optical counterpart, 3 1 have V 
magnitudes and V - I colors compatible with foreground main 
sequence stars. This fraction of X-ray sources (~ 3%) is com- 
parable to the estimated c ontamination in simil ar X-ray obser- 
vations (see, for example, iDamiani et alj|20 06). However, we 



will not exclude a priori these sources from our list of candidate 
cluster members, since they can be Class II YSOs with optical 
colors altered by effects due to their PMS nature (as disk gas 
accretion on stellar surface). Instead, they will be marked with 
a specific tag. 

2.2. Cross identification with GLIMPSE catalog 

GLIMPSE observations in the field containing NGC 6611 
cover the entire 33' x 34' WFI FOV, with the exception of a 
region at North- West of about 7' x 10'. Following the explana- 
tory manuafl we selected 41985 IRAC sources falling in this 
region. 



available at http://www.astro.wisc.edu/sirtf/docs.html 
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Table 1. Results of the cross-identification. 



Number of stars WFI detection 2MASS detection IRAC detection Number of X-ray sources 



146 no no no 146 

20476 no no yes 

2732 no yes no 26 

17768 no yes yes 64 

12535 yes no no 196 

1071 yes no yes 178 

2648 yes yes no 74 

2782 yes yes yes 370 



We identified the sources in common to both the GLIMPSE 
catalog and the catalog compiled by GPM07 with an identifi- 
cation radius of 0'3; this value is equal to the astrometric pre- 
cision of GLIMPSE catalog relative to 2MASS, as reported in 
the explanatory manual (in GPM07 we used the 2MASS Point 
Source Catalog as astrometric reference). In this work we also 
include 146 X-ray sources without any optical or 2MASS coun- 
terpart, which were not included in GPM07, to identify their 
possible IRAC counterpart. With the chosen identification ra- 
di us, we expect about 2 6 spurious identifications, evaluated as 



m 



Damiani et al. (2006). 



The results of the cross-identification are summarized in 
Table Q] We found 2782 stars with WFI, 2MASS and IRAC 
detections. This sample shows a significant spatial clumping 
in correspondence of the cluster and it includes 370 X-ray 
sources, suggesting that it is mostly composed by cluster mem- 
bers. 

Table Q] also shows that the 146 X-ray sources without 
any optical or 2MASS counterpart are not detected even in 
IRAC observations. Fig. Q] shows the spatial distribution of 
these sources, evidently clustered near the center of the ACIS 
FOV. They can be low-mass cluster members, with masses 
lower than the limit of our optical-IR data, since our catalog in- 
cludes X-ray sources down to the limiting magnitudes in WFI, 
2MASS and IRAC observations. However, it is also possible 
that some of these objects can be extragalactic sources, de- 
tectable in X-ray since the nebula is less dense in correspon- 
dence of the cavity cleared by the cluster itself. 

GPM07 used the detection in ACIS observation as mem- 
bership criterion. This was justified by the fact that almost all 
our X-ray sources with optical emission are PMS stars compat- 
ible with an age between 0.1 and 3 Myears. Fig. [2] shows the 

V vs. V - I diagram for the stars in WFI FOV with errors in 

V smaller than 0.1'" and in V - I smaller than 0.15'". In this 
diagram, diamonds mark the optical sources detected also at 
5.8fim. About 54% of these objects are bright stars (V < 16 m ), 
while the remaining 46% are fainter and mainly concentrated 
in the PMS region of the color-magnitude diagram. This re- 
gion overlaps with that traced by X-ray sources, confirming the 
youth of these optical-IRAC stars. The V vs. V — I diagrams 
of the sources detected in the other IRAC bands have similar 
characteristics. Hereafter, the color-magnitude and color-color 
diagrams of this paper will include only stars with errors in the 
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Fig. 1. WFI image, in the I band, of the central region of 
NGC 6611. Crosses identify the X-ray sources without both 
optical and infrared counterpart. The box outlines the 17' x 17' 
ACIS FOV. 



specific magnitude smaller than 0.1"' and in color smaller than 
0.15"'. 

In the next two sections we will describe the criteria 
adopted to select stars with circumstellar disk based on the 
excesses in IRAC colors, and we will compare them with the 
criterion based on 2MASS data, used in GPM07. These cri- 
teria are applied to stars with magnitude errors smaller than 
0.1 m ; this is a very stringent rule if applied to faint stars: at 
3.6/im more than 50% of stars fainter than 12.5'" have error 
greater than 0.1 m . The same 50% fraction is found at 12™ at 
4.5/jm and at 10.5"' at both 5.8/ira and 8.0/im. It is evident, 
then, that a criterion for the selection of stars with disk using 
all IRAC bands simultaneously does not allow to select faint 
sources with disk, since most of them have large errors in the 
IRAC bands at longer wavelenghts. We will partially overcome 
this problem using two independent disk diagnostics: the usual 
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Fig. 2. V vs. V - I diagram of stars in WFI FO V. The thick 
solid line is the ZAMS (from ISiess et al.1 12000b . at the dis- 
tance of 1750 pc and with Ay - 1.45'" (the average extinction 
of field stars, evaluated in GPM07). The dashed lines are the 
isochrones at 0. 1, 0.25, 1, 2.5, 3 and 5 Myears, with the average 
extinction appropriate for cluster members (Ay = 2.6 m ). The 
extinc tion vector is obtained from the law of Munari & Carrarol 
19961) . The diamonds are optical sources with emission at 5.8 



yum. 



color-color IRAC diagram (Sect. [3) and a diagnostic based on 
suitable reddening-free color indices (Sect. 0J, involving opti- 
cal+2MASS photometry and only a subset of IRAC bands. 

3. T-Tauri stars from IRAC color-color diagram 

The [3.6]-[4.5] vs. [5.8]-[8.0] diagram is an excellent tool 
to identify stars with circumstellar disk (see, for example, 
Allen et al.ll2004l) . In this diagram stars with photospheric col- 
ors are clustered around the origin, while stars with disk show 
colors significantly different from zero. The two populations 
can be easily distinguished even if the stars are affected by 
large interstellar extinctions. This is possible since the redden- 
ing vector is almost vertical (due to a very small reddening in 
[5.8]-[8.0]) or it points to bluer [5.8]-[8.0] colors (due to the 
partial overlap of the IRAC 8.0 fim band with the interstellar 
silicate feature), depending o n the adop ted extinction l a w; see , 
for example. iMegeath et all J2004 and iFlaherty et all J2007I) . 
Unlike what happens for other possible infrared diagrams, then, 
in the IRAC color-color plane the locus of reddened photo- 
spheres and that of sources with intrinsic red colors are not 
blended, even for large extinctions. 

With this diagram it is also possible to distinguish roughly 
embedded Class I and Class II T-Tauri YSOs (Young Stellar 
Objects): the former have both colors redder than the latter 
since the presence of the collapsi ng envelope surrou nding the 
star and the disk; see, for example, |AJlenet_al. (2004) who used 
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Fig. 3. Color-color diagram for IRAC sources in WFI FOV 
(dots). The box marks the approximate locus of Class II YSOs 
(crosses). The dashed lines roughly separate reddened photo- 
spheres (left), reddened Class II YSOs (center) and Class I 
YSOs (in the upper right part of the diagram, marked with 
triangles). Squares mark stars either classifiable as Class II or 
Class I YSOs. The r eddening vecto r s with Av=3 m and Ay=5 



are ob tained from Meg eath et al.1 (120041) . and IFlaherty et al 



d2007l) . respectively. 



the models developed b y b'Alessio et all dl998l Il999i l200lh 



for Cla ss II stars and by iKenvon et al. ( 1993b and Calvet et al 
d 19941) for embedded Class I stars. 



3.1. Classification of selected YSOs 

GLIMPSE data of the Eagle Nebula w ere already analyzed, i n 



combination with 2MASS data, by llndebetouw et al.l (2007). 



We reanalyze the data, using also different methods, for sev- 
eral reasons. First, we note that the selection of stars with K 



excess performed by llndebetouw et alJ (I2007I) and GPM07 is 
different. The latter used the reddening free color indices, a 
method much more efficient than the use of the T-Tauri locus 



(Mey er et all U997) in infrared color-color diagrams, adopted 



by the former. Moreover, in GPM07 only stars with color er- 
rors {(Tcoiors) smaller than 0.15'" have been used, and we ap- 
ply the same approach here in order to be more conservative 
as possible. With this conditi on we produce a d i fferen t IRAC 
color-color diagram respect to llndebetouw et aL ( 2007 ). In the 
present paper, this diagram is used together with reddening-free 
color indices defined with IRAC bands, in order to perform a 
selection of the cluster members with disk as much as possible 
model-independent. Furthermore, 2MASS and IRAC data are 
complemented by optical and X-ray data, that are fundamental 
to asses the nature of the candidate members. 
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Fig. [3] shows the color-color diagram for the IRAC sources 
in WFI FOV used to classify the YSOs. Stars with colors typ- 
ical of photospheres are clustered around the origin, with a 
spread consistent with photometric uncertainties and redden- 
ing, and are clearly separated from sources with intrinsic red 
colors. 

To select and classi fy stars with d i sk, we used the redden- 
ing vector obtained by Megeafh et al.l (120041) from the redden- 
ing law o f Mathis (1990), shown as the inclined vector in Fig. 
[3] The problem of the correct extinction law in IRAC bands is 
still open. For this reason, in Fig.[3]we also show the extinction 
vector (the Ak = 5'" arro w) obtained from the mean reddening 
law recently estimated by Flahert y et al.l d2007l) in the direction 
of five nearby star-forming regions. It is evident that the differ- 
ence between the two extinction laws is relevant only for a very 
large extinction. 

Using this diagram, we identify 147 Class II YSOs as the 
stars within the box in Fig. |3](crosses), taking into account pho- 
tometric uncertainties and extinction. This box was obtained 
from the models of stars with circ umstellar disk developed by 



D'Alessio et al. (1998, 1999, 2001). In addition, 13 stars are 



unambiguously classified as Class I YSOs (triangles in Fig . [3j . 

22 stars (marked with squares in Fig. O are sources with 
only one IRAC color redder than the Class II locus. Those with 
[3.6] -[4.5] > 0.8™ can be either Class I sources or very embed- 
ded Class II YSOs; those with [5.8] - [8.0] > 1.2 m are usually 
interpreted as Class II YSOs, with no emission detected from 
the inner region of the disk, due to a larger inner hole or to 
the disk inclinat ion with respect to th e fine of sight. However, 
as suggested by Kenvo n et al. ( 1993 ). also Class I YSOs with 
large centrifugal r adil and optically thin envelopes, that give 
a silicate band in emission, may have a spectral energy distri- 
bution compatible with the latter colors. For these reasons, all 
these 22 stars are not clearly classified. The possible contami- 
nation by extragalactic sources is discussed in the next section. 

The IRAC color-color diagram allows us to select a total 
of 182 members with excesses, with 120 new identifications 
with respect to GPM07, where BVIJHK bands were used. This 
result is not in disagreement with GPM07, since in that work it 
was not possible to classify most of these stars. Among these 
120 stars, in fact, only 48 have good measurements in K, 25 
in both V and / (shown in Fig. HI and just 9 in all these three 
bands simultaneously. We also note that 19 of these stars are 
also X-ray sources (among the total of 57 that are inside the 
ACIS FOV) and that the V - I and V values of the subsample 
with good WFI photometry are compatible with the cluster. 

Fig. [5] shows the IRAC colors of 76 X-ray sources with 
good photometry in the four IRAC bands. A large number of 
these sources (55) are classified as Class II YSOs while none 
as Class I sources. This is likely due to lower luminosity in 
X-ray of Class I objects respect to more evolved YSOs, in ac- 
cord to that showed by previous studies about star-forming re- 
gion s analyzed with X-ray observations deeper than our one 
(i.e. lPrisinzano et alJEooih . 
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Fig. 4. V vs. V-I diagram analogous to that in Fig. [2] Diamonds 
mark candidate stars with disk selected by the IRAC color- 
color diagrams but not in GPM07 
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2 the distance from the central star at which the infalling material 
from the envelope with most angular momentum, i.e. near the equato- 
rial plane, hits the circumstellar disk. 



Fig. 5. IRAC color-color diagrams of the X-ray sources 
(crosses). Dots mark the stars with disk selected from the IRAC 
diagram of Fig. [3] for the sky region falling in the ACIS FOV. 
The Class II locus and the reddening vectors are analogous to 
those in Fig. [3] 



3.2. Contaminating sources 

In order to evaluate the contaminati on due to extragalacti c 
sources, we use the criteria defined in lGutermuth et al.l (120081) . 
These criteria allow to identify candidate AGNs, or galaxies 
dominated by PAH emission, using various color-color and 
color-magnitude IRAC diagrams. Using these criteria, we iden- 
tify only 4 stars that can be extragalactic sources, but their Ay 
values predicted by the optical colors are too small (between 
2.5'" and 5'") to confirm this classification. 

However, in IRAC color-color plane PAH-rich (Polycyclic 
Aromatic Hydrocarbons) star-forming galaxies can be found in 
the locus at -0.1 m < [3.6] - [4.5] < 0.6" 1 and [5.8] - [8.0] > 
l m (in the part of the diagram redder than the Class II locus). 



6 



M. G. Guarcello et al.: Circumstellar disks in NGC 661 1 



O 

I 
I 



30:00.0 



35:Q0.0 



X^'X 



HP 



X 



-13:40:00.0 ' 







a 




- 45*0.0 

- 

- 50:00,0 



55:00.0 X 



■30.0 20.0 10.0 18:1 
-i i P 




■x^ 

X 



X 



30.0 



20.0 18:10.0 
. i . . . . I . . 



Fig. 6. IRAC image of the Eagle Nebula at 8.0 pm, centered on 
the cluster and covering approximately a region of 22' x 32'. 
The circles mark Class II YSOs, while the crosses mark Class I 
YSOs and stars compatible with both classifications. The small 
boxes mark massive stars. The box delimits the structures 
known as "elephant trunks". North is up, East is to the left. 

We have selected 25 candidate disk-bearing members in this 
region of the diagram. Some of them could be PAH-rich star- 
forming galaxies, and we will definitively classify them in our 
subsequent paper. 

3.3. Spatial distribution of Class II and Class I YSOs 

Fig. [6] shows the IRAC image of the Eagle Nebula at 8.0/im, 
with over plotted the Class II objects (circles) and the sources 
sorted either as Class I or as reddened Class II YSOs (crosses). 
In the center of the cluster (approximately in the center of the 
image) there is a lack of sources of the latter group, with respect 
to the outerregions. This is not a real result, since it is due to the 
intense diffuse emission of the dense structures in the central 
region of the nebula, like the elephant trunks inside the box in 
Fig. [6] that complicates the extraction of IRAC point sources. 
These structures, in fact, are heated directly by the radiation 
from the massive members of NGC 661 1, and are very bright 
in the IRAC bands at longer wavelengths. Besides, these diffuse 
structures are very dense, accounting for a very high extinction 
inside them. 



The evidence that star formation activity is still ongoing in 
the central region o f NGC 6611 was provided b y seve ral au- 
thors. For example, McCaughrean & Andersen! ( 2002h . using 
observations at high spatial resolution, showed that the trunks 
are stellar "nurseries", with several embedded YSOs inside 
them. The trunks are continuously eroded by the intense in- 
cident UV flux emitted by nearby massive cluster members. 
As results, some of the young sources formed inside the trunks 
are emerging from the photodissociation regions that delimit 
these structures. As effect of the intense diffuse emission by 
the nebula, we identify only one Class I YSO emerging from 
the trunks, at a = 18 : 18 : 52 and 6 = -13 : 49 : 38 (the cross 
inside the box in Fig. [6]). 

As shown in Fig.[6j 12 sources classified either as Class I or 
reddened Class II are clustered in a region at North- West. Only 
one of the few Class II stars present in this region has a faint op- 
tical counterpart (with V = 23.34'"). Therefore, this small clus- 
ter can be associated to a denser intra cluster medium and/or 
to more recent star formation events, if compared with the cen- 
ter of NGC 661 1. Th e presence of this rich sta r-formation site, 
already suggested by Indebetouw et al. ( 20071) . shows that the 
star formation activity is still ongoing in the whole nebula, and 
not only in the central region. 

4. Cluster members from reddening-free color 
indices 

In Sect. |2] we discussed the different sensitivities in the four 
IRAC bands and that this can limit the selection of stars with 
disk with the IRAC color-color diagram. To partially over- 
come this problem, we select the stars with infrared excesses 
also using four suitable reddening-free color indices, that will 
be described in det ail in the following (see also GPM07 and 
iDamiani et alj |2006 ) . 

4.1. Properties of the Q color indices 

Q indices are defined in order to compare infrared colors with 
V - I, the latter assumed to be representative of photospheric 
emission: 



Qvuyp] =(V-I)-(J- [sp]) x E V -i/Ej- 



[sp] 



(1) 



where Qvuyp] is the index; V, I and J are the magnitudes in 
the respective bands; [sp] is the magnitude in a specific IRAC 
band; Ey-i and Ej-[ sp ] are the reddening in the correspond- 
ing colors. EV_f ha s been obtained from the reddening law of 
Munari & Carraro (1996), while the extinction in IRAC bands 



has been taken from the reddening law o flMathisll 990( see Fig. 

The effectiveness of this approach is that Qvn[s P ] indices al- 
low us to determine excesses in stars for which the magnitudes 
in all the IRAC bands are not well measured simultaneously, if 
V, I and / magnitudes are known. 

The identification of stars with excesses is performed us- 
ing the diagrams Qvrjyp] vs - J ~ [sp\ in Figs. [7] In these di- 
agrams the extinction gives a shift only along the x - axis, 
while the excess in the [sp] band downward the y - axis. In 
this way, stars with excesses and reddened photospheres can be 
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Fig. 7. Diagrams 2hj[ !p ] vs. the J - [sp] colors for the stars in WFI FOV (points). Circles mark stars with excess in the corre- 
sponding index; the horizontal dotted-dashed lines are the lower limits forphotospheric indices; the inclined dashed lines separate 
the locus of normal stars from that of reddened sources (or stars with small excesses). Crosses mark the optical sources used to 
define these limits (as explained in the text). The squares mark stars with multiple identifications. 



separated, even if a significant number of stars remain unclassi- 
fied between reddened photospheres and stars with excesses. In 
Figs. [7] the stars with excesses are marked by circles; they are 
defined as the stars with the Qvnyp] indices significantly (i.e. 
by more than 3 times the error in Q indices) smaller than the 
horizontal dotted-dashed lines, that mark the lowest limits of 



photospheric Qvuyp] indices. We define these limits as the low 
boundaries of the loci, in these diagrams, of the IRAC sources 
with optical counterpart that are clustered around the origin in 
the IRAC color-color plane in Fig. [3] The spatial distribution 
and the positions in the other color-color and color-magnitude 
diagrams of these stars confirm that they are field normal stars 
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Table 2. The first column yields the Qvij[sp] indices used in this 
work; the second column shows the number of stars with excess 
detected with each Qvu[ S p] index; the third column shows the 
total number of stars for which we compute the indices. 



Q index 


Number of stars with excess 


Number of stars 


QviJ[3.6] 


113 


1746 


QviJ[4.5] 


124 


1474 


Gv/J[5.8] 


79 


357 


Gv/48.0] 


104 


297 



A total of 174 stars with excesses are selected 



or cluster members without disk. In Fig.[7jthey are marked with 
crosses, and it is evident that they are separated either from the 
stars with excesses and the reddened sources. Hereafter we will 
call NS locus that of the normal stars in the diagrams in Fig. [71 
EX locus that of the stars with excesses; UNC locus that of the 
stars compatible with both interpretations. 

The minimum excesses necessary so that a star with disk 
can be selected, by the use of the Qvij[sp] here and Qimass in 
GPM07, is mass and age dependent. To analyze these property 
of Q indices, we compute these minimum excesses in [3.6] and 
K for stars with masses equal to 0.5, 1, 1.5, 2, 2.5, 3, 4, 5 solar 
masses and with ages of 0.1, 0.5, 1, 2.5, 3.5, 5 Myears. Their 
photospheric indices are co mputed with the colors predicted 
by the evolutionary tracks of Siess et al. (2000), using the color 
transformation of Kenyon & Hartmannl ( 19951) . The minimum 
excesses are then the excesses in [3.6] or K necessary so that 
their Q indices become more negative than the lower limit of 
photospheric colors for more than 3cr, where cr is the mean 
error in Q of the sources in the WFI FOV. 

Fig.[8]shows the variation of the minimum excesses for one 
among the four indices used in GPM07 to detect excesses in K 
and for that used here to detect excesses in [3.6]. It is evident 
that for each mass and for each Q index the minimum excesses 
have irregular paths with increasing age. This is due to the ir- 
regular behavior of the indices when both the involved colors 
become more red (or more blue), as in comparing normal stars 
with different effective temperatures. 

Fig. [8] also shows that Qvijk is more efficient in detect- 
ing excesses in stars with lower masses respect to Qvij[3.6] in- 
dex. This behavior depends on the shape of the isochrones in 
the Q diagrams. For example, the minimum of the 5Myears 
isochrone in Qvu[3.6] of Fig. UJ occurs at ~B9 spectral type, 
corresponding to 2.5 M . This implies that stars of lower mass 
with small excesses easily fall in the UNC region, where we 
are unable to distinguish between small excesses and redden- 
ing and therefore we need stronger excesses to unambiguously 
asses the presence of disks. The corresponding minimum in 
Qvijk diagram occurs at ~K4 spectral type (using the 5 Myears 
isochrone), making this index much more efficient in detecting 
excesses in stars with smaller masses. 



4.2. Stars with Q excesses 

Table [2] gives the total number of stars for which we compute 
the indices and the number of stars with excesses. These num- 
bers can be lower limits since the criterion we adopt to detect 
the stars with excesses is very conservative. 

If we consider the total number of stars in our catalog with 
good photometry in the bands used to define the indices (third 
column in Table |2), we find a higher percentage of stars with 
excesses in the IRAC bands at longer wavelengths. This con- 
firms that a large number of optical sources are detected in the 
IRAC bands with longest wavelengths thanks to their nature of 
PMS stars, as suggested in Sect. 12.21 

Table |2] also shows a drop of the numbers of selected stars 
with excesses at [5.8], due to the smaller numbers of stars with 
good measurement in this band. This is an effect of the de- 
crease in sensitivity with increasing wavelength. Except for 2 
stars, the excesses in the three short IRAC wavelength bands 
are correlated. This suggests that they are related to the same 
physical region of the disk. This region can only be the inner 
rim at wh ich disk dusts sublimate. As propos ed by several au- 
thors (i.e. b'Alessio et alJll998i 1 1 999l 1200 lh . in fact, this re- 
gion is optically thick and it emits as a blackbody at the dust 
sublimation temperature (between 1500K and 2000K). 16 stars 
have excess only at [8.0]. This is a consequence of the fact that 
[8.0] can be affected by the 10pm silicate emission feature. A 
Class II YSO, then, can show excess in [8.0] even if the inner 
disk is mostly evacuated, as we can hypothesize for at least 4 
among the 16 members of NGC 6611 with excesses detected 
only with Qv«[8.0] (since the other 12 fall in the UNC loci in 
the other Qvu[sp] diagrams). 

5. Comparison among the used disk diagnostics 

Using the Qvu[sp] indices, we select 30 new stars with excesses 
(including 5 X-ray sources) not found with the Qimass indices 
in GPM07. Among these stars, as explained, 16 have excesses 
only at [8.0]; the other 14 have large errors in K band or fall in 
the UNC loci of the Qimass diagrams in GPM07. 

Similarly, 186 stars with excesses in 2MASS bands (mostly 
in K) selected in GPM07, are not classified as stars with ex- 
cesses with Qvu[sp] indices. This number of stars can seem 
very large, since usually more stars with excesses are observed 
at longer wavelengths. However, it must be noted that 94 of 
these stars have poor IRAC photometry, consistently with the 
different sensitivities between 2MASS and GLIMPSE catalogs, 
and therefore they are not included in IRAC sample. All the 
other 92 fall in the UNC loci of the diagrams in Fig. [7] 

We conclude that the set of indices defined in this work 
is consistent with that defined in GPM07 (with 144 sources 
with at least one excess in 2MASS bands and at least one in 
IRAC bands). This consistency means that if Q indices detect 
excesses in some particular 2MASS or IRAC band, the lack 
of excesses in the other bands is almost always due to poor 
photometry or to the ambiguity between excesses and redden- 
ing, as explained above. This is compatible with a scenario in 
which excesses in K (and H) are due to the same physical re- 
gion of the disk of the excesses in IRAC bands. It is important 
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Fig. 8. Trend of the minimum excesses necessary to identify stars with excesses in K and [3.6] for the studied masses and ages 
values. The excesses in K are detected with the use of Qvijk index, defined in GPM07; the excess in [3.6] with Qvu[3.6]- 
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in Fig. [7]; the stars with excesses detected in 2MASS but not in 
IRAC bands are marked as triangles. 



to note, however, that this consistency rules only between the 
Q indices, that are an efficient diagnostic for the selection of 
Class II YSOs having moderate disk inclination (with respect 
to the line of sight), because of the use of the optical and J 
bands. The color-color IRAC diagram is instead a peerless tool 
for the selection of embedded YSOs and Class II stars with 
highly inclined disk. Therefore the two methods may be con- 
sidered complementary. 

Fig.fTOland Fig.fTTIshow how the excesses detected with the 
IRAC color-color diagram translate into those detected with Q 
indices, and vice versa. In Fig. [TOl diamonds mark the 70 can- 
didate stars with disk selected in GPM07 with Qimass indices 
having good photometry in all IRAC bands, while dots mark 
the candidate stars with disk selected with the IRAC color-color 
diagram. Only 8 sources with excesses detected with Q2mass 
indices are close to the origin of the diagram (in the locus of 



Table 3. Number of stars with excesses detected with the dif- 
ferent diagnostics used in this paper and in GPM07. 



Disk diagnostics 


Numbers of stars with excesses 


IRAC color-color diagram 


182 


Qimass indices 


330 


Qvuup] indices 


174 



A total of 458 different candidate members with disk are selected. 



the normal stars), while the other have red IRAC colors. The 
characteristics of these 8 sources will be studied by SED anal- 
ysis in our subsequent paper. Fig. [TTI shows the Qvij[3.6] index 
of candidate cluster members with disk selected with the IRAC 
color-color diagram (marked with triangles). As expected, al- 
most no one of these stars fall in the NS locus. 

As further test, Fig.fHlshows the K - [3.6] vs. [3.6] - [4.5] 
diagram for the stars in ACIS FOV, with the X-ray sources and 
the stars with excesses in 2MASS bands selected in GPM07. 
All the stars with excesses in 2MASS bands have K - [3.6] > 
0.7 m and/or [3.6] - [4.5] > 0.2'", with the exception of a couple 
of stars that have also both IRAC colors ~ 0. 

In the diagram in Fig. [12] it is also clear that the distribu- 
tion of the stars with normal colors (marked with points) shows 
a gap at K - [3.6] ~ 0.5'", likely due to a rapid increase of 
the interstellar extinction at the distance of the Eagle Nebula: 
the foreground contaminating sources are clustered around the 
origin of the diagram and they are clearly separated from the 
sample of the more extincted stars, dominated by background 
sources. In fact, the gap between these two samples is popu- 
lated by a large number of X-ray sources, mostly associated 
with the nebula. 

Table [3] summarizes the number of stars with excesses se- 
lected with the diagnostics used here and in GPM07. 

6. Catalog of the candidate cluster members 

Combining all the diagnostics used in this paper and in 
GPM07, we select a total of 1264 candidate members of 
NGC 661 1: 790 candidate Class III members (inside the ACIS 
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Fig. 10. IRAC color-color diagrams of the stars with disk se- 
lected with the IRAC diagram of Fig.[3](points) and of the stars 
with excesses in 2MASS bands selected in GPM07 (diamonds). 
The Class II locus and the reddening vectors are analogous to 
those in Fig. [3] 



J-[3.6] 

Fig. 11. Diagram Qvij[3.6] vs. 7 - [3.6] for the stars in the WFI 
FOV, with overplotted the stars with excesses selected with the 
IRAC color-color diagram (triangles). 



FOV) and 474 candidate Class II and Class I PMS members 
(inside the WFI FOV). 

The catalog, available in electronic format, comprises: 

- stars IDs and coordinates; 

- the magnitudes and the errors in BVIJHK and IRAC bands; 
if some value is not available it is set equal to "NAN"; 

- a tag (tagx) that is equal to 1 if the star is also an X-ray 
source, otherwise the tag is equal to 0; 

- a tag (tagM) that is equal to A01 if the star is classified 
as a candidate members with disk only by its position in 
the IRAC color-color diagram, A 10 if the infrared excesses 
are detected only by Q indices and Al 1 if both diagnostics 



* + 



A",= 10 m 



0.2 0.4 0.6 

ra.aW4.5l 



Fig. 12. Color-color diagram with K and IRAC bands for the 
stars in ACIS FOV (points). Diamonds are stars with excesses 
in 2MASS bands detected with Qimass indices defined in 
GMP07, while crosses mark X-ray sources. The reddening vec- 
tor, corresponding to Ay - 10 m , has been obtained from the 
extinction law of lMegeath et alj ([2004). 



detect the emission from the disk; B if the star is a disk-less 
X-ray source; C if it is an X-ray source with optical colors 
compatible with a foreground main sequence star (see Sect. 
13. 



7. Spatial variation of disks frequency 

In order to confirm the results obtained in GPM07, we take ad- 
vantage of the new list of cluster members obtained here. Inside 
the central 17' x 17' ACIS FOV, where X-ray data are available, 
we select 790 candidate PMS members without disk, thanks to 
their X-ray emission and the presence of an optical/IR counter- 
part, among which 31 are likely foreground contaminants. In 
addition, in the same FOV we found 257 stars with NIR ex- 
cesses, 54 more than in GPM07, that are candidate to have a 
circumstellar disk (1 18 are also X-ray sources). 

The average disk frequency inside the ACIS FOV is 24% + 
2%, larger th a n the value obtained in GPM07 (19% + 1%). 
Olivei ra et a j] (120051) found a disk frequency equal to ~ 58% 
in a small region inside the ACIS FOV. The different fraction 
are likely due to the different sensitivities of the tw o surveys. 
In fact, the L survey used in lOliveira et al.l (120051) is deeper 
than the GLIMPSE catalog ( about 1 magnitude res pect to [3.6] 
sources catalog). Moreover, lOliveira et alJ (12005b shown that 
the disk frequency increases with decreasing mass of the cen- 
tral star, likely since the more massive is the central star the 
more rapid is the erosion of the circumstellar disk, as it was 
found in other star for ming regions (for example, the study of 
Carpenter et al. , 2006, on the Scorpius O B association). This 
may be the origin of the discrepancy with lOliveira et a l. (2005) 
results. 
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Fig. 13. Percentage of members with disk in the ACIS FOV vs. 
the estimated incident flux from OB members. The four values 
are 31% + 4%, 27% + 4%, 21% + 3% and finally 16% ± 3%, 
from left to the right. 



However, more than the absolute value of disk frequency, 
we are interested to its spatial variation respect to the position 
of massive members, that allows to understand if UV radiation 
emitted by the latter members may alter the disk lifetimes of 
nearby T-Tauri YSOs. In order to reach this goal we need only 
to be sure to use a consistent criterion in different sky positions. 
GPM07 already used this approach, and we verify here their 
results using the new list of cluster members. 

For every cluster member, with and without disk, we com- 
pute the incident flux emitted by 52 massi ve members of the 



cluster, with spectral class earlier than B5 dHillenbrand et al 



19931) . To this calculation, we use the projected distances from 
massive stars, as already explained (Sect. [TJ. We then compute 
the disk frequencies for various bins of incident flux, as shown 
in Fig. [13] It is evident that the main GPM07 result, i.e. that 
members with disk are more frequent at larger distances from 
massive stars, where they are irradiated by lower UV fluxes, is 
confirmed. Disk frequency, in fact, increases with decreasing 
incident fluxes: from 31% ± 4% in the bin with lowest flux, to 
27% + 4%, 21% + 3% and finally 16% + 3% in the bin with 
highest flux. In the histogram in Fig. [13] the bins size are de- 
fined in order to have the same number of disk less members 
in each bin. These disk frequencies are obviously more reliable 
than those computed in GPM07, thanks to the new and more 
complete list of cluster members. This paper, then, reinforces 
the result of GPM07 about the influence of massive stars on the 
evolution of nearby members with disk in NGC 661 1. 

7.1. Spatial distribution of the cluster 

The panels in Fig.[14]allow to compare directly the spatial dis- 
tribution of the cluster members and the intensity of the radia- 
tion emitted by massive members. The cluster members are not 
equally distributed in the ACIS FOV, with an empty region at 
South- West. The intensity map shown in the upper right panel 
of Fig. [14] is more concentrated that the members density, re- 
flecting the central concentration of massive members. 



The two lower panels of Fig. [14] show the spatial distri- 
butions of disk-bearing and disk-less members. Note that the 
latter stars have a more symmetric distribution; while the for- 
mer are more present where the UV flux is low, and also their 
distribution present an "hole" in the center of the field, where 
the UV flux is lower. Moreover, the region corresponding to 
the highest fluxes of Fig. Q~3] have a radius of 0.6 parse cs, the 
distance from massive stars at which lBalog et al. I d2007h found 
a drop in disk frequency in NGC 2244. 

The upper right panel of Fig. [T41 a lso allows a more easy 
confrontation of our result with that of IBalog et alj (12007b . As 
explained, these authors studied the spatial variation of the disk 
frequency respect to the positions of massive stars in the young 
cluster NGC 2244. They have not found a correlation between 
the two distributions, but observed a significant drop in disk 
frequency for distance smaller than 0.5 pc from O stars. The 
central regions of the intensity map in Fig. [14] corresponding 
to the highest fluxes in Fig. Qj] has a radius equal to about 0.6 
parsec, so the decrease of the disk freque ncy in the la s t bin i n 
Fig.[13]is in agreement with the finding of lBalog et al. ( 2007 ). 

To have an estimation of the extension of NGC 66 1 1 , even 
if it has not a symmetric spatial distribution, we use the 2- 
parameters density profile of King ( 19661) : 



aif) 



0~o 



1 + (r/r core y 



(2) 



where <tq is the central cluster density, r is the distance of 
stars from cluster center and r core is the core radius. We cal- 
culated from the observed radial density profile: ctq = 149 
+8W starspc an d r core = 1.39 ± 0.08pc, in agreement with the 
results obtained in GPM07. 



8. Summary 

In this paper, we analyze Spitzer/IRAC data of the young open 
cluster NGC 661 1, in the Eagle Nebula, in a large field wide ap- 
proximately 30' x 30'. In a previous paper, we have already se- 
lected the members of this cluster using our published B VIJHK 
and X-ray multi band catalog. The obtained list of cluster mem- 
bers, both with and without disk, is used to verify that the 
former are more frequent at larger distances from the massive 
members of the cluster. 

In this work, in order to select new members with disk, 
we use and compare two different disk diagnostics: the IRAC 
color-color diagram, that use all the IRAC bands simultane- 
ously, and four suitable reddening-free color indices (Q in- 
dices), each defined to detect the excesses in a specific IRAC 
band. 

Using the IRAC color-color diagram, we identify 182 
Young Stellar Objects (147 Class II; 13 Class I and 22 with 
colors compatible with both classifications). This diagnostic is 
very efficient to select largely embedded young sources, but 
it can be used only for stars with good photometry in all IRAC 
bands. We discuss how significantly this affects the selection of 
faint YSOs. With the use of Q indices we partially overcome 
this problem, since they allow us to select stars with excesses 
in each IRAC band taken alone; with this diagnostic, we select 
174 YSOs with excesses in IRAC bands. Among them, 66 stars 
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Fig. 14. The upper left panel shows the spatial distribution of candidate cluster members. The dotted square delimits the ACIS 
FOV. The upper right panel shows the intensity map of the radiation emitted by the massive members of the cluster. The different 
shades correspond to the four bins in Fig. [13] The lower left panel shows the spatial distribution of candidate members with disk; 
the lower right panel that of candidate members without disk. In all the panels, the cross marks the center of the ACIF FOV 



cannot be detected with the IRAC color-color plane, since they 
are not well measured in all the IRAC bands. 

Combining the outcome of both diagnostics, we identify 
146 new cluster members with disk with respect to our previ- 
ous work, which was based on 2MASS photometry alone (116 
from the IRAC color-color diagram and 30 from Q indices). At 
the end, our catalog includes a total of 474 candidate Pre-Main 
Sequence members with disk and 790 without disk (the latter 
only in the central 17' x 17' field, where we found 118 X-ray 
sources with disk). 

Comparing all the Q indices defined with 2MASS and 
IRAC bands (with the exception of [8.0] band), we claim that 
they are all sensitive to the emission from the same physical re- 
gion of the disk, namely the inner rim at the dusts sublimation 
radius. These indices will be efficiently used, then, to select 
Class II stars that are not largely embedded (i.e. for which it 
is possible to detect the photospheric emission) and whose in- 
ner disk is not evacuated. The embedded YSOs associated to 
the cluster can be selected only thank to the IRAC color-color 
diagram, and not with the Q indices. 

We discuss evidence suggesting that star formation activity 
is ongoing in the outer region of the Eagle Nebula, and not only 
in the center of NGC 6611 as showed by previous works. For 



example, the diagnostic based on IRAC color-color diagram 
allow us to identify a probable new star formation site, rich of 
Class I and embedded Class II objects, at North with respect to 
the center of the cluster. 

In the central 17' x 17' field we find an average disk fraction 
equal to 24% + 2%. In this field, the disk frequency also has a 
strong dependence on the incident radiation emitted by mas- 
sive members, varying from 31%+ 4% to 16% + 3% across the 
whole range of values of incident flux. This is an evidence of 
the influence of massive stars radiation on the evolution of cir- 
cumstellar disks and star formation process, as already shown 
in GPM07. 
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